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The electrical conductance of atomic metal contacts represents a powerful tool to detect nano- 
magnetism. Conductance reflects magnetism through anomalies at zero bias 1-7 - generally 
with Fano lineshapes - due to the Kondo screening of the magnetic impurity bridging the 
contact. 8 ' 9 A full atomic-level understanding of this nutshell many-body system is of the 
greatest importance, especially in view of our increasing need to control nanocurrents by 
means of magnetism. Disappointingly, zero bias conductance anomalies are not presently 
calculable from atomistic scratch. In this Letter we demonstrate a working route connecting 
approximately but quantitatively density functional theory (DFT) and numerical renormal- 
ization group (NRG) approaches and leading to a first-principles conductance calculation 
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for a nanocontact, exemplified by a Ni impurity in a Au nanowire. A Fano-like conductance 
lineshape is obtained microscopically, and shown to be controlled by the impurity s-level po- 
sition. We also find a relationship between conductance anomaly and geometry, and uncover 
the possibility of opposite antiferromagnetic and ferromagnetic Kondo screening - the lat- 
ter exhibiting a totally different and unexplored zero bias anomaly. The present matching 
method between DFT and NRG should permit the quantitative understanding and explo- 
ration of this larger variety of Kondo phenomena at more general magnetic nanocontacts. 

The electrical conductance of metals through magnetic nano-contacts, including STM tip 
tunnelling 16 and mechanical break junctions (BJ) 7 displays, in analogy to quantum dots 10 , a zero- 
bias anomaly reflecting a Kondo many body effect. 8 ' 9 First principles calculations of the contact 
conductance and of the anomalies are currently unavailable. The standard DFT based Landauer 
method 11 is mean-field in character, and thus invalid in presence of many body effects. Alterna- 
tively, the NRG 12 solution of Anderson impurity models (AIM) 13, 14 is known to yield correctly the 
zero-bias conductance 15 : but the AIMs are treated as barely more than toy models, their parameters 
generally guessed or fixed phenomenologically, rather than based on first principles. 

Our approach starts from a first-principles DFT calculation but, instead of ending as usual 
with the standard calculation of the mean-field conductance 16 ' 17 , it continues with the building of 
geometry-dependent AIMs, each explicitly including the channels, orbitals, and especially sym- 
metries suggested by the DFT electronic structure. The key point is that the numerical parameters 
are fine tuned by requiring these AIMs to yield in mean field the same phase shifts suffered in 
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Figure 1 : Electronic structure of a Au wire with a Ni impurity in bridge and substitutional geome- 
tries, a) DFT projected density of states (PDOS) at the Ni s orbital in bridge geometry. Blue dotted 
line: spin up; full black line: spin down, b) PDOS at the Ni d zx orbital in bridge geometry. Black 
dashed line: LDA+U down spin PDOS with U = 3eV; c) electronic bands of a clean Au chain; 
d) PDOS at the Ni d zx , d zy orbitals in substitutional geometry. The effect of U is negligible in this 
case. 
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DFT by a spin polarized traveling electron of each symmetry. We exemplify directly our method 
by straightforward application to a prototype Au-Ni-Au monatomic chain contact, a test system 
chosen because of its simplicity: in fact, Au possesses just one s-like conduction band, and the 
Ni bridging impurity is free of complications involving anisotropy and spin-orbit coupling, which 
affect other transition metal atoms. 4 

We start off from the recent DFT calculation of the nanocontact electronic structure 17 and 
examine the symmetry- selected projected density of electron states (PDOS). The Ni atom in the Au 
chain (axis parallel to z) has two low energy geometries: bridge (B) and substitutional (SUB) (left 
and right insets in Fig. [I]). In B (preferred at zero stress) the PDOS indicates a well defined Ni 3d zx 
spin down empty orbital (Fig jTJ?). The Au chain conduction channel is a single 6s nondegenerate 
band crossing the Fermi energy E F (Fig. [T£), hybridized with both Ni 4s and Ni 3d zx orbitals - 
the first even under reflection across the Ni xy plane, the other odd. Electrons traveling in the Au 
"leads" will thus scatter differently off the Ni impurity in the two symmetry channels, even (e) and 
odd (o). The DFT calculated magnetic moment M = 2/j,b{S z ) is 1.30/xb, enhanced over the bare 
value of a S = 1/2 impurity state by co-polarization of nearby Au atoms, but dropping in fact to 
~ 1.0/iB after correction of DFT through a Hubbard-f/ term (so-called LDA+U) which raises the 
3d zx energy 17 (dashed curve in FigjTJ?). Summing up, the B geometry implies electron hopping 
across Ni 3d zx which has odd symmetry and S ~ 1/2, as well as across Ni 4s which has even 
symmetry and is very little spin polarized. In SUB geometry (of higher total energy 17 , but probably 
accessible at large stress), the electron states can be labeled by m, the orbital angular momentum 
about the z axis. The empty spin-down Ni state orbitals are now a 3(d zx , d zy ) degenerate pair with 
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\m\ = 1. Orthogonality to the m = Au conduction band results in a PDOS 5-like peak just 
0.02 eV above E F (Fig. |T]d). This suggests 5 = 1; and the DFT calculated Ni magnetic moment 
is in fact M = 1.91/ie, rising to 2.02/i£ upon addition of a U term. 17 Summing up, in the SUB 
geometry the Ni impurity has S ~ 1 and the magnetic orbitals are "spectators". The Au conduction 
electrons hop onto the Ni 4s orbital in the even symmetry channel, whereas they find no Ni valence 
orbital in the odd channel. 

The DFT calculations also provide the electron transmission matrix through the impurity, 
yielding directly the Landauer ballistic conductance. However, in standard DFT the moment is 
artificially frozen, breaking spin rotational invariance, thus wrongly predicting a different conduc- 
tance for spin up and spin down electrons 17 - in zero external field there simply is no telling what 
is up and what is down. Through the Kondo effect, 8 symmetry is restored by screening of the spin 
by the conduction electrons. This many body effect is usually studied by applying a viable spin- 
symmetric technique such as the NRG 12 to a microscopic, but phenomenological AIM. This wastes 
the detailed information provided by DFT, unless some way is found to preserve it. The extraction 
of Anderson model parameters from DFT has a long history. 18 Very recent work on Fe impurities 
in bulk Au and Ag also addressed this problem. 19 Our proposed joining step between DFT and 
NRG is the construction of a different, specifically designed Anderson model for each symme- 
try and geometry, with parameters embodying quantitatively the DFT calculated electron-impurity 
scattering amplitudes and phase shifts. For each symmetry (e or o) an incoming electron wave with 
spin initially along x will, after scattering on Ni (magnetic moment parallel to the z-axis), rotate in 
the x — y plane, clockwise if the coupling is FM or counter-clockwise if it is AFM. Within DFT, 
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Angle (rad) 


Bridge DFT (HF-AIM) 


Subst. DFT (HF-AIM) 




-0.28 (-0.28) 


-0.19 (-0.19) 


60 


1.10(1.17) 


-0.12 (-0.00) 



Table 1: Deflection angles. Spin rotation angles e / o (at the Fermi energy) of an incoming 
transverse-polarized conduction electron wave after scattering at the Ni impurity for even and 
odd conduction channels. DFT ab initio angles are compared with the AIM Hartree Fock ones 
(in brackets), optimized by adjusting the AIM parameters. Positive e / o correspond to antiferro- 
magnetic electron-impurity coupling, negative to ferromagnetic. o substitutional is zero in the 
AIM model, as there is no Ni empty valence orbital of this symmetry, and inter-site exchange is 
neglected. 

electron-impurity scattering leads to calculable spin rotation angles, e / o = 2(8K — 5L ), which 
measure the phase shifts and thus the coupling sign and strength of the Ni impurity spin to the Au 
conduction electrons. The DFT calculated rotation angles (Table [T]) show that in geometry B the 
coupling is large and AFM in the odd channel, small and FM in the even channel. In geometry 
SUB the coupling is on the contrary weakly FM in both channels. 

We now have all the ingredients needed to construct the minimal AEVIs which describe the 
Au-Ni-Au contact. The models must include the Ni 4s orbital, with energy e s and creation operator 
s^ , and the 3d a orbitals, with energies e a and creation operators S aa , where a denotes spin and 
a = xy, xz, . . . labels orbital symmetry. Electron transport in the Au leads takes place in the 6s 
band (the 5d band is filled) leading to Au conduction state combinations with p = e, o spin-cr 

6 




Figure 2: Conductance properties in the bridge geometry, a) NRG calculated Kondo temperature 
T K illustratively shown as a function of the on-site d zx interaction U . The physical value for Ni 
3d zx is indicated by the arrow, b) Color plot of the zero bias conductance as a function of temper- 
ature (normalized to the Kondo temperature) and of the effective s orbital energy e s /D treated as a 
parameter for illustrative purposes. The physical value for Ni 4s is indicated by the horizontal line. 
When e s >> Ep the low temperature (T < Tr) conductance is zero, and increases by moving e s 
towards E F , hitting a maximum at e s ~ —D, declining below -D. c) Low temperature differential 
conductance as a function of bias voltage V s d applied to the junction. The color of curves refer to 
the value of e s in the left bottom panel. Note the Fano lineshapes, due to variation of the even-odd 
channel interference with varying e s . (See supplementary material). 
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creation operators c J kpa and energy e p k- The AIM Hamiltonians have the general form 

H = E e pk c{ pa c kpcT + Vks (cL* s * + H- c ) 

kp<r kcr 

+ E E *W { C lp* d aa + H.c) (1) 
kpcr as 

+ E C * S Ua + E 6 « 4<A„ + flint [S, d a ] . 
a act 

In |T|), Vks and Vk pa are hybridization parameters between leads and impurity orbitals, while 
Hint [s, d a ] includes the mutual interactions between impurity orbitals. All parameters in ([TJ) de- 
pend on the geometry of the contact, and must reflect the physics provided by the DFT calculation. 
Our practical procedure is to estimate first approximate AIM parameters from direct inspection of 
the DFT electronic structure PDOS, and to operate successively a finer adjustment by comparison 
of the Hartree-Fock calculated AIM spin rotation angles with those directly provided by DFT. Be- 
cause the AEVIs include only a few orbitals and interactions, the agreement is semi-quantitative: 
yet it can, as we will show, be made very significant. 

Bridge - In B geometry, the Ni 3d zx orbital is directly hybridized with the odd conduction 
states, the Ni 4s with the even states, and all other Ni 3d states can be dropped. The interaction H int 
includes a Hubbard U ~ 2eV (deduced by the exchange splitting of spin minority and majority 
3d zx PDOS peaks ), and an intra-atomic s-d FM exchange Jh ~ —0.3eV (extracted by the spin 
splitting of Ni 4s, co-polarized by the magnetic 3d xz ). A band-impurity hybridization width T zx = 
2npoV x 2 z ~ 0.25 eV (po ~ O.leV^ 1 is the conduction band density of states at Ep) is deduced from 
the broadening of the DFT calculated PDOS Ni 3d zx peaks (see FigjTJ?). Moreover, the majority 
and minority peak asymmetry around E F is strongly reduced by addition of a realistic Hubbard U 
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in LDA+U (see dashed peak in FigjTJ?), so we may assume a particle-hole symmetric e zx = —U/2, 
with U ~ 2 — 3 eV. Finally, Ni 4s is a broader orbital, for which U ~ and T s ~ 3 eV, as extracted 
from the PDOS linewidth. The s-level energy parameter e s is at the outset totally uncertain. The 
fine adjustment to reproduce the DFT spin rotation angle yields e s = —0.4 eV, just below E F . 
For the sake of illustration, since e s will sensitively control the Fano lineshape of the conductance 
anomaly, we choose nevertheless to explore different values of e s , both above and below Ep, as 
might be relevant in other nanocontacts. Note that the bridge configuration has an additional direct 
Au-Au hopping bypassing the impurity (see inset of FigTiJ)), its effect taken equal to T s (for further 
detail see supplementary). 

The AIM Hamiltonian ([TJ with the parameters adjusted as above yields in the HF approx- 
imation spin rotation angles quite close to those calculated by fully realistic DFT (Table [T] see 
also supplementary material), so that the low energy description of the B configuration is close 
to being quantitatively correct. The next step is the NRG treatment of model ([T]) (see Methods). 
At low temperatures, Jh becomes irrelevant and the e and o channels behave independently. The 
conductance G = G sin(5 e — 5 ) 2 , where G = 2e 2 /h, is calculated by following, during the 
NRG flow, the phase shifts 5 e / Q (see Methods). The competition between the Kondo effect in the 
strongly AFM odd channel and the resonance in the weakly FM even channel results in the zero 
bias conductance plotted in Fig. [2J5 as a function of temperature and of e s , the conventional posi- 
tion of Ni 4s relative to E F . At large unphysical e s >> D (D = 2.5eV is half bandwidth), both 
channels acquire a n/2 phase shift, with a destructive interference and zero conductance. At the 
fitted value of — D <C e s < 0, the conductance is instead expected to be G/G ~ 0.5. This deviates 
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considerably from the DFT calculated Landauer conductance, relatively close to the unitary limit 
in both spin channels, G = G up + Gd OW n ~ 0.5 + 0.4 ~ 0.9G - 17 The large difference under- 
scores, as anticipated, the impact of Kondo correlations, ignored by DFT. By further decreasing 
e s < 0, the conductance reaches its maximum value G/G = 1 and then drops. This effect of e s 
can be assimilated to that of a gate voltage, controlling the conductance lineshape versus source- 
drain voltage. Following Ref. 15 we calculated the nonequilibrium conductance as a function of 
the source-drain voltage V sd (details given in supplementary material). The conductance shows 
a Fano-like resonance near zero bias, as shown in Fig. |2p, and as expected on general grounds. 9 
The large predicted Kondo temperature ~ 100 K, see Fig. [2^, should make this anomaly readily 
accessible in a break-junction transport experiment. 

Substitutional - In SUB geometry, magnetism occurs in 3d zx and 3d zy Ni orbitals that do 
not hybridize with the Au 6s conduction band. The interaction term H int [s, d a ] in Eq. ([TJ) should 
now include all multiplet exchange splittings between 3d zx and 3d zy degenerate magnetic orbitals. 
The intra-atomic ferromagnetic exchange with the 4s Ni orbital is weak, see the small DFT rotation 
angles in Tab{TJ and all the electron-impurity phase shifts are now FM. Unlike regular AFM Kondo, 
the ferromagnetic Kondo exchange is well known to be irrelevant at low temperatures, 8 where 
the magnetic orbitals asymptotically decouple from the conduction band. Only Au 6s and Ni 
4s hybridization will effectively survive, leading to a tt/2 phase shift in the even channel. The 
zero bias conductance will thus approach unitarity G(V = 0) < G , with a broad Breit-Wigner 
lineshape of width T s expected for the spectral function. While confirming that, the explicit NRG 
calculation with the Anderson model (1) with orbitals and parameters appropriate to the SUB 
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geometry shows in addition that the FM Kondo antiscreening mechanism gives rise at the Fermi 
level to a logarithmically sharp spectral function"pimple", see Fig.[3](right panel) and also Ref. 20, 
pushing G/Gq up to unity at zero bias. This is again different from the DFT Landauer mean 
field conductance, 17 G T + Gj. ~ (0.40 + 0.38)G ~ 0.8G . The lineshape difference of this 
expected "ferro Kondo" anomaly in SUB geometry and of the Fano lineshape of ordinary Kondo 
in B geometry, respectively right and left panels in Fig. |3j is striking. More generally, no ferro 
Kondo zero bias anomaly appears to have been reported so far anywhere, so its existence in the 
present or in other nanocontacts is open for experimental detection. While more work will be 
needed to describe the effect of a magnetic field in a ferro Kondo case, the effect of temperature 
should be to cut off the logarithmic cusp, modifying the conductance in the form 1 — const/ In 2 T. 
The final conceptual aspect of spin rotational symmetry is worth mentioning. As is well known, 21 
in the regular AFM Kondo, for example of our B geometry, the impurity approximate 5 = 1/2 
is absorbed by Kondo screening into an exact singlet, a trivially rotationally symmetric state. In 
the FM Kondo coupling of SUB geometry the impurity's approximate 5 = 1 is turned by Kondo 
antiscreening into an exact triplet, fully SU(2) rotationally symmetric, endowed with its 25 + 1 
degeneracy. 

Summarizing, this work represents a first attempt at joining together the quantitative advan- 
tages of DFT electronic structure with the correct NRG many body Kondo physics of a magnetic 
impurity in a break junction metal contact. Because the intermediate Anderson models can only 
represent approximately the full electronic processes, the procedure is semi-quantitative. It is 
nonetheless revealing, showing in the regular Kondo case an ab initio example of Fano lineshape; 
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Figure 3: Conductance vs. source-drain voltage of regular and of Ferro Kondo in the two geome- 
tries. Sketch of the different zero bias anomalies expected for the current through an Au-Ni-Au 
junction in the B (left panel), with realistic e s = —0.4 eV, and SUB (right panel) geometry. Oppo- 
site to the Fano lineshape calculated in the B geometry (see also Figj2]:), the Kondo antiscreening 
in SUB geometry yields a spectral function (see inset) which, for small energies, exhibits a loga- 
rithmic "pimple" (see text). 
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involving phase shifts and Kondo processes that are sometimes ferromagnetic and involve antis- 
creening; and predicting very different zero bias anomalies depending on impurity geometry and 
spin. The method should be applicable more generally, whenever well defined intermediate AIMs 
a) can be built and b) are soluble. On the experimental side, numerous magnetic atom nanocontacts 
have been realized using STM 3 6 . Mechanical break-junction studies have just begun, so far ex- 
ploring a different type of phenomenon, namely the emergence of Kondo effects in ferromagnetic 
metals 7 , but many more should be possible in the future, when the kind of antiferro-ferro Kondo 
switching described here should hopefully be pursued. 
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METHODS 



DFT We carried out DFT calculations within the spin-polarized generalized-gradient approxima- 
tion (a-GGA) with the Perdew-Burke-Ernzerhof parametrization 22 for the exchange and corre- 
lation energy, with same parameters as in Ref.17. All the calculations were done with the plane 
wave PWscf code, included in the QUANTUM-Espresso package. 23 Ultrasoft pseudopoten- 
tials 24 were employed and kinetic energy cutoffs of 30 Ry and 300 Ry were used for the wave 
functions and the charge density, respectively. For both B and SUB geometries, the Au-Au bond 
length was taken 2.80 Ato avoid spurious magnetizations due to DFT self-interaction errors, 17 
whereas all the Au-Ni distances were fully optimized. The calculation supercell consisted of 16 
Au atoms plus one Ni atom periodically repeated in all three directions. The wire-wire distance in 
the xy plane (perpendicular to the wire axis z) was 10.58 A, making spurious interactions between 
periodic replicas negligible. In the B configuration, the Ni atom lies in the x z plane. Convergence 
with respect to A;-points and smearing parameters was carefully checked. 

Transmission and reflection amplitudes, and the spin rotation angles (see the Supplementary 
material) were calculated using the Choi and Ihm's method 11 generalized to ultrasoft pseudopo- 
tentials, 25 as implemented in the PWCOND code (a part of the QUANTUM-Espresso package). 
The self-consistent potential in the first part of the supercell described above, of length equal to 
the Au-Au distance, was used to build the periodic potential of the left and right leads, while the 
potential in the rest of the supercell was used as the scattering region. The zero bias ballistic con- 
ductance was obtained using the Landauer-Biittiker formula 26 from the transmission coefficient at 
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Ep (with all spin moments frozen). Spin-orbit effects were not taken into account in the present 
study. The PDOS in Fig. [Tjwere calculated directly from the scattering states. 

NRG In our NRG code (for a review see Ref. 12) we implement the U{\) charge symmetry (quan- 
tum number Q, the total charge with respect to one electron per site and orbital) and the SU (2) 
spin symmetry (quantum number S, the total spin). We choose the Wilson discretization parameter 
A = 2 (or 1.8 when calculating the Green functions), and keep up to 1500 states per iteration when 
calculating dynamical quantities, or 800 if we are only interested in the energy spectrum. The 
Kondo temperature can be expressed as 8 T K = where w = 0.4128 is the Wilson coefficient, 
T the hybridization linewidth in the magnetic channel, and Z the quasiparticle residue. The lat- 
ter can be extracted from the self-energy according to: Z^ 1 = 1 — (dYi(iu n ) / d iu n )\ iuJn=Q . The 
self-energy is £ = — G^ 1 + Gq 1 where G and G are the impurity Green functions calculated by 
NRG, respectively in the presence and absence of interaction. The resulting Kondo temperature in 
our model is shown in Fig. [2] as a function of U. In the bridge geometry the zero bias conductance 
is evaluated as a function of temperature by direct inspection of the NRG flow. It can be expressed 
as G = 2e 2 sin 2 (5 e — 5 Q )/h. The phase shift 5 e / Q are related to the two lowest energies of states 
with quantum numbers (Q, S) = (1, 1/2), which correspond to the cost of adding an even or odd 
electron to the ground state that has quantum numbers (0, 0). We calculated the difference between 
these phase shifts, hence the zero-bias conductance, as a function of the temperature T extracted 
from the NRG iterations. We note however that, while we are quite confident about the values at 
low temperatures, those at high temperatures must be taken with caution, since the spectrum is still 
far from a Fermi liquid one. In the substitutional geometry we calculate the Ni 4s spectral function 
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following the guidelines of Ref.27. 
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